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ABSTRACT 



1. INTRODUCTION 
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The detection of 1.809 MeV gamma-ray line emission 
from the Cygnus X complex by the COMPTEL tele- 
scope is one of the most convincing proves of massive 
star nucleosynthesis in our Galaxy. The Cygnus X 
complex is an extremely active nearby region, con- 
taining several star forming regions, OB associations 
and young open star clusters. It houses some of the 
most massive stars known in our Galaxy and concen- 
trates a large number of Wolf-Rayet stars. Thermal 
radio continuum emission and intense [C n] line emis- 
sion reveals widespread ionisation, and at least 60 
individual H n regions have been identified. In order 
to understand the 1.809 MeV gamma-ray line emis- 
sion from the Cygnus X complex, and to compare 
the observations to theoretical nucleosynthesis calcu- 
lations, we modelled the multi-wavelength spectrum 
of the region by means of an evolutionary synthesis 
model. Our investigation leads us to the following 
conclusions: 

1. Stellar wind ejection is the dominant mechanism 
for the observed 26 A1 enrichment in the Cygnus 
X region 

2. Cyg OB2 is by far the dominating massive star 
association in Cygnus X and 1.809 MeV line 
emission from 26 A1 produced in this association 
should be detectable by the spectrometer SPI 
on INTEGRAL 

3. There is only low supernova activity in the 
Cygnus X complex and the 60 Fe lines should be 
below INTEGRAL 's detection sensitivity 

The detectability of an individual massive star clus- 
ter (Cyg OB2) by INTEGRAL provides a fantastic 
opportunity for future nucleosynthesis studies using 
gamma-ray line spectroscopy. We will explore the 
scientific potential of such observations, and provide 
estimates for gamma-ray line intensity distributions 
based on the massive star census of the Cygnus X 
region. 
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Cygnus X designates one of the most prominent 
extended emission regions in galactic radio contin- 
uum surveys, located between I = 70° — 90° and 
b = ±5°. At high angular resolution, the e mission 



breaks int o hundreds of continuum sources (Wend- 
ker, 1984), most optically invisible or heavily red- 



dened, largely owing to foreground obscuration by 
the Cygnus Rift at a distance of 700 pc. 

1.809 MeV gamma-ray line emission fr om the Cygnus 



X reg ion has been first reported by Del Rio et al 



(1996) based on the analysis of 2 years of COMP- 



TEL observations. They observed a diffuse extended 
emission feature with a total 1.809 MeV line flux of 
7 10 -5 p h cm^s' 1 . To under stand the origin of this 
feature, |Del Rio et al. (1996| ) modelled 26 A1 nucle- 
osynthesis by adding the expected individual con- 
tributions from known Wolf-Rayet (WR) stars and 
supernova remnants (SNRs) in this area. In that 
way they predict a total 1.809 MeV flux of 4.2 10~ 5 
ph cm _2 s _1 from which ~ 85% is attributed to WR 
stars. A close inspection of their list of candidate 
sources reveals, however, that this flux estimate is 
dominated by a few nearby (< 1 kpc) objects. Since 
distance estimates for individual WR stars and SNRs 
are quite uncertain, it is questionable how accurate 
such a candidate source census can be. 

In this work we present an alternative approach to 
model 26 A1 production in Cygnus which is based on a 
census of massive star populations in the field (such 
as OB associations and young open clusters). The 
basic parameters, like distance and age, are deter- 
mined much more reliable for stellar ensembles than 
for individual objects, hence we believe that our 
model is more accurate in the prediction of mas- 
sive star nucleosynthesis. On the other hand, mas- 
sive stars that are not situated in associations are 
not accounted for by our model, so strictly speak- 
ing we only obtain a lower limit. Globally, how ever, 
most massiv e stars are located in associations ( |Gar- 



many, 1994 ), hence the underestimation should at 
most amount to 20% (we will show later that it prob- 
ably is even smaller) . Additionally, our approach will 
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be much l ess sensitive to comp leteness limits than the 
model of Del Rio ct ah (1996 ), since the population 
parameters can be estimated reasonably well without 
knowing all members of the population. In particu- 
lar, we also estimate nucleosynthetic contributions 
from stars that already faded away by means of an 
evolutionary synthesis code (Cervino et al., 2000; see 
also Cerviho et al., these proceedings). Although 



yield (in units of Mq) using 



Del Rio et al. (1996) also try to estimate this con- 
tribution from known radio supernova remnants, it 
is clear that this approach probably misses an im- 
portant fraction of recent 26 A1 production due to in- 
completeness of SNR catalogues. In particular old 
remnants are difficult t o detect aga inst the galactic 
background radiation ( |Green, 1991 ), but their 26 A1 
still contributes to the today 1.809 MeV radiation. 



2. FREE-FREE EMISSION 



The radio continuum emission from the Cygnus X 
region is primarily optically thin thermal radiation 
that arises fr om individual H n regions and diffuse 
ionised gas (Wendker, 1970). Ionised gas is also 
nicely traced by the far-infrared transition of [C n] 
at 158 /Ltm, and the correlation of diffuse [C n] and 
thermal radio emission in Cygnus X is further ev- 
idence for t he presence of a bundant ionised gas in 
this region (Doi ct al., 1997). 



Thermal radio emission has also shown to be impor- 
tant in the context of 26 A1. Knodlseder et al. (1999 ) 
have demonstrated that the 1.809 MeV intensity dis- 
tribution follows closely the galactic distribution of 
free-free emission, mapped for exampl e by the DMR 



instru ment aboard COBE at 53 GHz ( Bennett et al 



1992). In fact, this correlation is one of the strongest 
pieces of evidence that 26 A1 production arises from 
massive stars which also power the free-free emission 
by ionisation of the interstellar medium. To express 



26 A1 p roduction in terms of ionisation, Knodlseder 



|(1999| ) proposed to normalise the 1.809 MeV emis- 
sion on the number of ionising photons produced by 
an representative ionising star (e.g. an 07V star), 
leading to an average equivalent 07V star yield of 
Y 26 7 V = (!-0 ± 0.5) 10~ 4 M for the entire Galaxy. 

From the COMPTEL 1.809 MeV 7 years all-sky map, 
Pliischke et al. (200C| ) derived a total gamma-ray line 
flux of $i. 8 09 Mcv = (7.9 ± 2.4) 10~ 5 ph cm^s- 1 for 
the Cygnus region, in rea sonable agreement wi th the 
earlier determination by Del Rio et al. (1996) from 
only 2 years of data. Integrating the intensity in the 
DMR 53 GHz free-free emission map over the same 
region results in a total flux of $53 GHz = 5400 ± 1200 
Jy, where the quoted error reflects the uncertainty 
that arises from the background level determina- 
tion and the subtraction of an underlying galactic 
ridge emission that is probably not associated to the 
Cygnus X complex. The ratio between both flux 
measurements is then converted into the equivalent 



Y° 7V = 7.9110 3 x 



1.809 MeV 



(ph cm 2 s 1 ) 



53 GHz 



(Jy) 



(i) 



Thus, for the Cygnus X region we obtain = 
(1.2 ± 0.4) 10 -4 Mq, compatible with the galactic 
average value. 



3. MASSIVE STAR POPULATIONS 



Our massive star census of the Cygnus region is b ased 
on the WEBDA catalogue of |Mcrmilliod (1998|) for 
open clusters and the compilations of Glarmany 



Btcnccl (1992) and Humphreys (1978) for galac- 
tic OH associations. The data have been comple- 
mented by information from the SIMBAD database, 
r ecent analyses of stella r associations by Hipparcos 
(Dc Zeeuw et al., 1999), and a recent i nvestigation 



of th e stellar population of Cyg OB2 (Knodlseder 
2000| ). For the investigated field 60° < / < 110° and 
\b\ < 15°, our final database contains 10 associations 
and 21 young open clusters, where young means an 
age of less than 50 Myrs, corresponding to the life- 
time of a 8 Mq star (i.e. the lowest initial mass that 
is believed to lead to a type II supernova). We find 
a total of 182 O stars within the OB associations, 
while the clusters only contain 19 O stars. 

To estimate the fraction of massive stars in Cygnus 
that is included in our association census, we ex- 
tracted from the SIMBAD database all stars in the 
field that are classified as type O. Following; the 
discussion of Garmany et al. (1982 ) this sample 



should be complete to ~ 2 — 3 kpc, corresponding to 
the maximum distances of massive star associations 
found towards Cygnus X (see Fig. [l]). An exception 
might be the heavily reddened Cy g OB2 association 
that contains about 120 O stars ( |Knodlseder, 2000 ) 
that are only partially identified in SIMBAD. Cor- 
rccting for this underestimation, we find a total of 
223 classified O stars in the field from which 185 (or 
83%) are found in one of the associations or clus- 
ters of our database. The remaining field O stars 
are mainly found towards Cep OB2 and north-east 
of Cyg OB3 showing a wide spread in distances be- 
tween 1 — 5 kpc. Thus we believe that our database 
present a fairly complete description of the massive 
star census in Cygnus, and that the small (17%) frac- 
tion of field stars and their large distances will not 
lead to a severe underestimation of the nucleosyn- 
thetic production in this area. 



For each of the associations or clusters in our 
database we estimated the distance by the method of 
spectroscopic and photometric parallaxes. H-R dia- 
grams were built from spectroscopic and photometric 
data from which cluster ages have been determined 
by isochrone fitting. Particular care has been at- 
tributed to the estimation of age and distance uncer- 
tainties since they directly tra nslate into uncertain - 
ties of our model predictions (Cervino et al., 2000). 
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Figure 1. Distribution of OB association (circles) 
and young open clusters (diamonds) in the Cepheus- 
Cygnus region projected onto the galactic plane (the 
positive x-axis points towards the galactic centre). 
The grey-shaded area delineates the boundaries of the 
investigated field, the labelled arc indicates galactic 
longitudes. 



We also derived initial mass functions from the data 
to estimate the initial richness of the populations (a 
detailed description of the methods will be presented 
in a forthcoming paper). 

To illustrate the spatial distribution of the investi- 
gated populations, we used our distance estimates to 
project the associations and clusters onto the galac- 
tic plane (see Fig. [l]). The most prominent feature 
in this representation is an elongated structure run- 
ning from the solar vicinity towards I ~ 70° to a 
distance of ~ 3 kpc, known as the local spiral arm. 
From the 201 O stars in our database 200 lie within 
associations and clusters of the local spiral arm, il- 
lustrating the extremely young age of this structure. 
60% of the O stars lie within a single association, 
Cyg OB2, situated at a distance of ~ 1.6 kpc. The 
extreme richness of this association, as revealed by 
near-infrared obser vations, has recently led one of us 
(Knodlseder, 2000) to identify Cyg OB2 as a young 
globular galactic cluster. From isochronc fitting we 
estimate the age of this cluster to 2.5 ± 1.0 Myrs. 



4. MODEL RESULTS 



Using our evolutionary synthesis model (Cerviho et 
al., 2000; see also Cerviiio et al., these proceedings), 
we estimated the integrated gamma-ray line and free- 
free radiation fluxes for the entire Cygnus region. In 
addition to 26 A1 we also included 60 Fe nucleosynthe- 
sis in our model to predict also the gamma-ray flux 



Table 1. Model predictions versus observations. 
Gamma-ray fluxes are given in units of 10~ D 
ph cmT 2 s~ x . 

predicted observed 



$1,809 MoV 


1.8- 


5.4 


7.9 ±2.4 


$53 GHz (Jy) 


2900- 


9300 


5400 ± 1200 


YO e 7V (1O- 4 M ) 


0.2- 


1.0 


1.2 ±0.4 


$1,137 McV 


0.1 - 


0.4 





Initial mass function slopes between T = — 1.1 and 
— 1.35 have been explored which covers the range 
of values that we found for the associations in our 
database. Our model predictions are summarised in 
Table [j]. For each quantity of interest we predict a 
range of possible values that includes uncertainties 
in distance, age, richness, evolution, and IMF slope 
for each association. We deliberately did not in- 
clude uncertainties related to theoretical nucleosyn- 
thesis yields, since we consider our model as a hy- 
pothesis that we want to test against observations of 
gamma-ray line emission. In total we predict a 1.809 
MeV gamma-ray line flux between (1.8 — 5.4) 10 -5 
ph cm" 2 s _1 where the lower value corresponds to 
the lower flux limit for an IMF slope of T = —1.35 
while the upper value is the upper flux limit obtained 
for a flat slope of T = —1.1. On average, about 80% 
of the 1.809 MeV emission in our model originates 
from 26 A1 ejected by stellar winds while only 20% 
comes from core-collapse supernovae. Also, OB as- 
sociations play a dominant role since they contribute 
more than 90% of the 1.809 MeV emission in Cygnus. 

The most prominent source of 26 A1 is probably Cyg 
OB2 for which we predict a median 1.809 MeV flux of 
8.5 10 -6 ph cm _2 s _1 , at the detection limit of cur- 
rent COMPTEL observations and above the sensi- 
tivity limit of SPI for a two-weeks observation. The 
flux uncertainty for Cyg OB2 is quite large, ranging 
from 1.1 10~ 6 to 2.3 10~ 5 ph cm' 2 s _1 , which mainly 
reflects our adopted age uncertainty of 2.5±1.0 Myrs. 
At this age the most massive stars of the association 
start to evolve off the main sequence which gener- 
ally means a considerable increase of the mass-loss 
rate, either in a LBV or a WR phase, leading to an 
enrichment of the interstellar medium with freshly 
produced 26 A1. Consequently, a small variation in 
the adopted age at this stage leads to a quite im- 
portant variation in the number of mass-losing stars 
(LBV and WR), which translates directly in a large 
variation in the interstellar 26 A1 enrichment. The 
existence of a few WR stars in Cyg OB2 indicates 
that, if formed coevally, the association should not 
be younger than ~ 2 Myrs. On the other hand, 
Cyg OB2 can not be much older than 3-4 Myrs since 
the absence of supernova remnants within the clus- 
t er boundary indicate s that all stars are still alive 
(Wendker et al., 1991). However, star formation in 



Cyg OB2 was prob ably not strictly coeval ( Massey 



from the radioactive decay of this longlived isotope. 



Thompson, 1991), hence the indicated age uncer- 
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Figure 2. Predicted median 1.809 MeV intensity 
distribution. 



tainty merely reflects the observed age-spread in this 
cluster. In this case, our median flux value corre- 
sponds to the age averaged 1.809 MeV intensity. 

To illustrate the spatial distribution of 1.809 MeV 
emission predicted by our model we generated an in- 
tensity map based on the median flux values that we 
distributed within the association boundaries using 
Gaussian shaped density profiles (cf. Fig. |2|). This 
figure may be compared to Fig. 1 of Pliischke et 
al. (these proceedings) which shows the COMPTEL 
1.809 MeV maximum entropy image of the Cygnus X 
region. The emission maximum in the COMPTEL 
map lies only 3° from the centre of the Cyg OB2 clus- 
ter, which is certainly not a proof, but a strong indi- 
cation that Cyg OB2 is indeed a prominent source of 
26 A1 (note that the angular resolution of COMPTEL 
at 1.809 MeV amounts to ~ 4° FWHM). Additional 
contributions in this area come from the OB associ- 
ations Cyg OBI, Cyg OB3 and Cyg OB9 which are 
situated west of Cyg OB2. 

However, the predicted 1.809 MeV flux from Cygnus 
is only marginally compatible with the observed flux, 
a problem that has already been e ncountered in the 
earlier work of Del Rio ct al. (199(f ) . We possibly un- 
derestimated the massive star population in Cygnus 
X since we did not included stars lying outside associ- 
ations, yet our discussion in Sect, ^suggests that our 
census is fairly complete. In particular, our model 
reproduces the observed 53 GHz free-free emission 
reasonably well, and an increase in the number of 
massive stars would quickly lead to an overproduc- 
tion of ionising photons in the Cygnus region. Ad- 



ditionally, the predicted equivalent yield Y^ 7 v falls 
also below the observed value, indicating that the 
problem is more likely related to an underestima- 
tion of the 26 A1 yield per massive star. The origin 
of the low Y26 7V values could lie in erroneous age 
determinations, since Yf^ 7 v depends strongly on the 
evolutionary state of the associations. Indeed, some 
of the OB associations (such as Vul OBI, Cyg OB7, 
Cyg OB8, Cyg OB9 and Cep OBI) show only poorly 
defined main sequences which makes it difficult to 
estimate their actual age. It is even ques tionable if 



the se associ ations are physically related (De Zeeuw 



et al., 1999) which questions the validity of treat- 



ing their 26 A1 production by means of an evolution- 
ary synthesis code. However, the number of O stars 
situated in these associations is fairly small (about 
17% of the total number), and it is unlikely that 
they alone are responsible for the discrepancy. Fi- 
nally, there remains the possibility that current nu- 
cleosynthesis models underestimate 26 A1 production 
in massive mass-losing stars. Indeed, a comparison of 
the observed galactic 26 A1 mass with a simple model 
for galactic 26 A1 enrichment also indicates a possible 
underestim ation of the globa l galactic 26 A1 yield by 
30 — 50% (Knodlseder, 1999). Recent calculations 
of 26 A1 nucleosynthesis in rotating massive stars sug- 
gest that ro tation could co nsiderably enhance 26 A1 
production ( Meynet, 1999|) , providing an interest- 
ing mechanism that could bring the theory in better 
agreement with the observations. 

The predicted 1.137 MeV flux arising from the de- 
cay of 60 Fe is rather small, and a detection of this 
emission by SPI seems highly improbable. The low 
flux is a result of the low supernova activity in the 
Cygnus X region, which is the dominant production 
channel for this radioisotope. On the other hand, the 
observation of large 1.809 MeV fluxes in conjunction 
with low supernova activity is a remarkable situation. 
Although we have no observational evidence that al- 
lows us to generalise such a correlation for the entire 
Galaxy, this indicates that supernovae may play a 
less important role for galactic 26 A1 production than 
previously thought. 
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